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ABSTRACT 

The first-generation stars in the ACDM universe are considered to have formed in dark halos with 
total masses in the range ~ 10 5 — 10 7 M Q at z ~ 20 — 50. These stars expected to be very massive and 
in some cases they end their lives as the first supernovae (SNe). We explore the problem of whether 
star formation in low mass dark halos 10 7 M©) was triggered or suppressed by the SN feedback 
from neighboring star-forming halos. We take into consideration mainly two effects by the SN shock: 
one is the evacuation of gas components from the halos and the other is the promotion of H2 formation 
because of the enhanced ionization degree by shock heating. Combining above effects, we find that the 
star formation activities in the neighboring dark matter halos (M ^ 1O 7 M0) are basically suppressed 
in case they are located close to the SN center, because of the gas evacuation effect. The critical 
distance within which the gas is blown away falls within the range ~ 0.3 — 1.5kpc depending on the 
SN energy and the halo mass. In addition, we find there is very little window in the parameter space 
where star formation activities in dark halos are induced or promoted by neighboring SN. 
Subject headings: cosmology: theory — hydrodynamics — stars: supernovae: general — ISM: super- 
nova remnants 



1. INTRODUCTION 

Star formation in the early universe has played a crit- 
ical role in subsequent evolution of the universe. First- 
generation stars may have substantially contributed to 
the cosmic reionization and metal pollution of the uni- 
verse by their radiative/kinetic feedback effects. These 
events are important for the formation and evolution 
of protogalaxies. The studies of the initial collapse of 
primordial pre-galactic objects in the ACDM universe 
have been done in b oth analytical (|Nishi fc Susa 1 119991: 
Tegmark et al.lll997l) and numeri c al method (lAbel et al.l 
19981 iFuller fc Couchmaf] 120001 lYoshida et a!T T2003). 



These objects are formed with masses of the order of 
~ 10 5 -10 7 M Q at redshifts z ~ 20-50. Theoretical stud- 
ies suggested that in these low-mass halos, H2 molecules 
are formed up to the level of ~ 10 -4 , because the virial 
temperature and the central density of t hese halos are 
high enough to activate H2 formation (jNishi fc Susa I 
1999) so that the cooling time becomes short enough. 
This small fraction of H2 is sufficient to cool the gas, 
which leads to the formation of first-generation stars in 
these halos. 

Early studies on the formation of primordial stars 
have been done in almost 



lYoshida et al.l 120061 ). These studies consistently 
suggested that first-generation stars are very massive 
(~ 30 — 500M Q ). Because of their extreme mass scale, 
first-generation stars emit copious amount of ionizing 
radiation, as well as a strong flux of ^-dissociating 
Lyman-Werner (LW) band radiation. Therefore, the 
radiation from the first stars dramatically influences 
their surroundings, heating and ionizing the gas within 
a few kiloparsecs around the progenitor star. 

If the primeval star is very massive (> 100M Q ), the 
ionization front from the star will break out of par- 
ent halos up to 10 7 M© in mass. Half of the baryons 
in the halo will be swept up into a dense shell that 
grows to the virial radius of the halo (~ 100 pc) by 
the end of the life of the star. Therefore, the gas den- 
sity within the shell radius will drastically decrease to 
low uniform densities o f 0.1 — lcm~ 3 prior to the su- 
pernova (SN) explosion (lAbel et al J 120071: lAlvarez et al 
2006t iKitavama et alj 12004 IWhalen. Abel fc Normal 
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)96l) . On the other hand, after 1990s, there have 
been a number of numerical studies of the formation of 
prim o rdial st ars (|Abel et alj 
19991 l200T 
1999: 
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The dark halos are also affected by the radiation 
from the first star in the neighboring halos. The lo- 
cal radiative feedback effect is sensitive to the distance 
from the source star, as well as the mass and the evo- 
lutionary stage of the target halo. The ionizing ra- 
diation from the star totally or partially photoevap- 
orates them before the SN blast ever reaches them 
(lAhn fc Shapirdl2007l:lMesinger et aT]|20fjllO'shea et all 
120051: ISusa fc Umemural 120061: IWhalen et all l2008aft . 
Therefore, these halos will be less massive or be gone 
altogether by the time the SN shock strikes them. It is 
also should be noted that the LW band radiation from 
the source star photodissociates H2 molecules in these 
halos i n case the halos are located close enough to the 
source (Gl over fc Brand ll200ll ISusall2007T) . 
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Some of ma ssive primordial sta r s end their lives as en- 
ergetic SNe. iHeger fe Wooslevl (|2002f ) suggested that 
the progenitor star of 10 - 40M Q dies in Type II SN 
blasts while those of 140 — 260M© dies in pair in- 
stability SN. On the other hand, it is important to 
note that some POP III stars lying in between the 
ranges indicated above may also explod e as hyper- 
novae (jTominaga. Umeda fc Nomotol l2007h . The stud- 
ies of first SN explosions have been performed in both 
SPH (jBromm. Yoshida fe Hernquistl 120031: iGreif et all 
20071) and grid co des (jKitavama fe Yoshidal 120051 : 
Whalen et alTl2008bh . In less massive halos (M < 
10 7 Mq), the ejecta first expands into a nearly uniform 
rarefied ionized medium, then interacts violently with 
the dense shell swept up in the progenitor H u region. 
Part of the energy of the blast wave is reflected into the 
center as a reverse shock and the rest pushes the shell for - 
ward (iKitavama & Yoshidal [2001 1 Whalen et al.ll2008bD . 
In massive parent halos (> 10 7 Mq), however, the SN 
remnant (SNR) will only expand 20 — 40 p c into the halo 
and then recollapse (j Whalen et al.|[2b08bf ). never reach- 
ing any nearby halos because the /-front does not break 
out these parent halos prior to the SN explosion. 

The shock wave from an SN that break away from 
the host halo impact and go through the neighboring 
halos those survived photoevaporation process. The 
shock wave may blow off the gas component; however, 
it also promote the formation of H2 and HD molecules 
in the gas because of the enhanced electron abundance 
by th e shock heating (jFerraral 119981: iJohnson fc Bromml 
2006t iKang fc Shapirol Il992t iNishi fc Susa I 119991: 
Oh fc Haimanl l2002t IShapiro fc Ka ngj 119871: IS usa et al l 
19981 : lUehara fc Inutsukall2000| ). We also note that the 
relic H 11 region including the photoevaporated halos 
recombines out of equilibrium, which is basically same 
physical condition of the postshock enhancement of 
H2/HD formation. The SN shock can sweep up this 
H2/HD and carry it into the nearby halo. Therefore, 
these molecules in the fossil H 11 region are important for 
the formation of secondary stars as well as those would 
be formed in the SN shock by collisional ionization. 

Star formation could be triggered by SNRs in several 
ways. Primordial SNRs may sweep up a dense shell of 
ambient gas or collide with the dense shell swept up by 
an H 11 region, either of which becomes contaminated by 
metals in the ejecta, subsequently comes to be dominated 
by its self-gravity. As a result, the shell is expected to 
fragment into smaller filaments/co res where populations 
of less massive stars are formed ( Machi da et al.l 120051: 
Mackev et aHl2003l : lSalvaterra et alJl2004t IWhalen et al l 
20081J 1. On the other hand, in rather massive halos 



(> 1O 7 M ), the H 11 region generated by the progeni- 
tor star is confined well inside the virial radius and the 
gas is kept neutral. In this case, first SNRs expand in 
neutral halos, heavily mix their interiors with heavy el- 
ements, and then rec ollapse without escaping the halo 
(jWhalen et al.ll2008bfl . 

The interaction of the first SN with neighboring halos 
in the early u niverse has not been investigated in detail so 
far. Recently. I Greif et al.l (|2007f ) performed a cosmologi- 
cal simulation of first SN explosion. They found that the 
shock from the SN can accelerate the star formation pro- 
cess in neighboring rather massive halos in which stars 
could be formed without feedback effects. They also indi- 



cated that if the SN explosion occurs in an H 11 region, the 
SNR comes to pressure equilibrium at half of the radius 
of the relic H 11 region because of its relatively large tem- 
per at urej^JTJiisJimits^ 

1.5kp c (jKitavama et all 120041; IWhalen Abel fc Normanl 
I2004h . In addition, ICen fc Riquelmd (|2008f ) performed 
numerical simulation of SNRs ram-pressure stripping 
cosmological halos. They focus on how heavy elements 
from the remnant are mixed with the halo gas. 

However, the number of halos studied by numerical 
simulations is restricted. More analytical criteria are 
needed in order to obtain a systematic understanding of 
first SN feedback effects on neighboring dark matter ha- 
los. In this paper, we investigate feedback effects by first 
SN onto the nearby halos in the early universe, which is 
potentially important for the total star formation activ- 
ities in the early universe. We use analytic arguments in 
order to obtain the universal criteria for feedback effects 
by first SN on the neighboring dark halos. The outline 
of this paper is as follows. The initial setup and the de- 
scription of SNR evolution are given in Section 2. In 
Section 3 and 4, the SN feedback effects are listed and 
quantified. Section 5 and 6 are devoted to discussion and 
summary. 

2. MODEL 

We consider SN explosions in the first collapsed ha- 
los (< 10 6 M Q ). The gravitational potential of these 
halos are so shallow that the ionizing radiation from 
the progenitor f irst stars can sweep out the gas of 
the host halos (lAbel et alJ 120071: lAlvarez et al| l20Q(i: 



Kitavama et~afl |2004 IWhale n. Abel fc Norma 



Wise fc Abell 120071) . As a result, the subsequent SNe 



can easily break away from the halo because of the de- 
creas ed gas density by photoioniz a tion prior to the explo - 
sion (jKitavama fc Yoshidal 120051 : IWhalen et afll2008bf) . 
In order to understand the nature of the SN shock ex- 
panding into an essentially uniform and ionized inter- 
galactic space, we assume spherical symmetry and ini- 
tially uniform ambient gas density of averaged cosmolog- 
ical density, po- Strictly speaking, this assumption is not 
valid since the gas density of the halo still slightly higher 
than po even after the feedback of UV radiation. This 
assumption should have some effects to increase the en- 
ergy of the SN shock, but we use this for simplicity. The 
energy loss mechanism from the SNR is dominated by 
bremsstrahlung for a first few years, then line emission, 
and then inverse Compton s c attering becomes imp ortant 
(jKitavama fc Yoshidal [20051 : IWhalen et~aTI l2008bl) . We 
take into account all of the cooling rate stated above 
as well as H2 cooling at low temperature. We use the 
fitting formula for these rat es fro m the compila t ion b y 
iFukugita fc Kawasakil (fl99l and iGalli fc Pallal (fl998ll . 
respectively. 

We consider the neighboring dark matter halos with 
total mass of 10 5 Mq ^ Mdh ^ 1O 7 M , whose baryonic 
fraction in mass is assumed to be the cosmic mean value, 
^b/^M- We have another free parameter in our calcula- 
tions, which is the distance from the SN center to the tar- 
get halo. Additionally, we use two typical fixed redshift 
when the SN explode (zsn = 20) and when the nearby 
halos virialized (z V ir = 30). We tested this issue with 
three different SN explosion energy E$n — 10 51 ,10 52 , 
and 10 53 erg. Throughout the paper, we work with the 
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ACDM universe with 17m 
fl h = 0.05. 



0.3, fl A = 0.7, h = 0.7, and 



2.1. Timescales 

We introduce three timescales i s , i CO oi, and tg, that 
characterize the important physical processes. They rep- 
resent the sound crossing timescale of the gas compo- 
nent in the neighboring halo, the cooling timescale, and 
the free-fall timescale, respectively. The sound crossing 
timescale of the gas component in the neighboring halo 
is described as 

*.(0 = A (i) 

where I denotes the length that we are interested and c s 
is the sound speed of the gas, respectively. This timescale 
is also interpreted as the expansion timescale of the shock 
heated gas if we substitute the size of the shock-heated 
region for I and use the sound speed of shocked gas as c s . 
The cooling timescale is defined as 



tcoolCF) 



ikT 



( 7 -1)A(T>,/h 2 ) : 
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where T and n are the temperature and the number den- 
sity, k and 7 denote the Boltzmann constant and the 
ratio of specific heats, respectively. /h 2 represents the 
fraction of H 2 , and A (erg cm~ 3 s" 1 ) denotes the cool- 
ing rate of gas. The free-fall timescale is written as 
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Here, G is the gravitational constant and p- 
virial density given by p v ir = 187r 2 p cr , where p, 
W~ 29 h 2 



(1 



^vir - 

Zvh) 3 g cm" 3 . 



vir is the 
1.9 x 



2.2. Evolution of SNR 

The time evolution of an SNR in intergalactic medium 
(IGM) is mainly described by following four stages: 

(1) The free- expansion stage. The free-expansion stage 

lasts until the SN ejecta sweeps up roughly the 
same amount of mass as their own in the sur- 
rounding medium. In this stage, the velocity 
of the SN ejecta decre ases linearly with radius 
dTruelove fc Mckedll999l) . 

(2) Sedov-Taylor adiabatic expansion stage. The expan- 

sion of the shock front is well approximated by the 
Sedov-Taylor solution. 

(3) Pressure-driven expansion stage. After the post- 

shock gas is cooled via the radiative cooling, the 
SNR depart from an adiabatic expansion. Geomet- 
rically thin shell is formed just behind the shock 
front. The expansion of the shocked shell is driven 
by the high pressure of the hot cavity. 

(4) Momentum-driven expansion stage. In this stage, 

the shocked shell expands conserving its momen- 
tum. 

In our calculation, we do not take into consideration the 
stage (1) because the duration of this stage is very short 



and it hardly affects the entire result. We also do not 
consider the stage (4) in IGM since the shock velocity at 
this stage is too small to activate H2 molecule formation 
or to evacuate the gas from the halo. We consider the 
scenario, in which the SNR gradually sweeps up mass 
and then collides with the neighbor halo. In fact, the 
remnant first collides violently with the dense shell swept 
up in the progenitor H 11 region, but we leave this effect 
for future works, since introducing such effect makes the 
analysis complicated. 

2.2.1. The Sedov-Taylor adiabatic stage in IGM 

The expansion of the shock front in the Sedov-Taylor 
adiabatic stage is describ ed by the self-similar solution 
(|Sedovill946l:lTavlorill950D . The radius and the expansion 
velocity of the shock front is written as 



R s = 1.15 



SN 



po 



1/5 



.2/5 



dR s 
dt 



0.460 



V Po 



1/5 



t 



-3/5 



(4) 
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where -Bsn and po represent the SN explosion energy 
and the ambient density, respectively. The postshock 
temperature is derived from the Sedov-Taylor solution 
and the Rankine-Hugoniot relation as follows: 



Tps : 



: 0.423 ( 



k ) 
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(6) 



where mu and p denote the atomic mass unit and the 
mean molecular weight, respectively. These solutions are 
based upon the assumption, that the gas is adiabatic, 
however, the postshock gas forms a dense shell which is 
cooled by the radiative cooling subsequently. After the 
cooling timescale of the shell becomes shorter than the 
expansion timescale of the shell (t CO oi < Rs/vs), the SNR 
move on to the pressure-driven expansion stage. 

2.2.2. The pressure-driven expansion stage in IGM 

The shock front expands by the high pressure in the 
hot cavity. We assume that the pressure inside the cavity 
decreases adiabatically as 



Pr, 



Pi 



Rs 
Ri 



-37 



(7) 



where P ca denotes the pressure inside the cavity, and also 
Ri and Pi represent the radius of the shock front and 
the postshock pressure at the beginning of the pressure- 
driven expansion stage. Combining this equation with 
the equation of motion of the shell, we have the sho ck 
radius/ velocity as follows (Sakashit a fc Ikeuchi Hl99"6T ). 
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We also obtain the postshock temperature as follows: 




2.3. Metal cooling 

In the present paper, we do not take into account the 
cooling rate due to heavy elements. In fact, the shock- 
heated gas in the nearby halos is expected to be pol- 
luted by metals ejected from the SN, although the abun- 
dance of heavy elements in the shock-heated gas is hard 
to evaluate. The ultra high resolution simulations would 
be indispensable to assess the degree of metal mixing. 
Aside from such difficulty, we can evaluate the metal 
abundance assuming complete mixing between the SN 
ejecta and the surrounding material. Based upon such 
assumption, the mean metallicity of t he swept up mass 
by the SN shock is Z /Z P) ~ 1CT 2 - 5 (|Greif et alJ [20071 : 
ISalvaterra et al.ll2004l ). In addition, the metallicity of 
the high-z IGM observed by Lya absor ption systems i s 
at a level of Z/Z Q ~ 10~ 3 to 10" 2 (|Songailal l200l . 
Thus, Z/Z & ~ 10 -3 to 10~ 2 could be a rough stan- 
dard to assess the effects of metal cooling. At such 
low metallicity, th e radiative cooling rate is h ardly af- 
fected above 10 4 K (Bochrin ger fc He nslcr 198j|), since it 
is dominated by H and He cooling. The cooling rate 
below 10 4 K is basically propor t ional to the amount 
of m e tals dBeniamin et alJ [20011: iDalgarno fc McCrav I 
19721 iRavmond et al l 1197a IRavmond fc Smith! 119771: 
Sutherland fc Dopitalll993l ). however, again at such low 
Z, the metal cooling rate is c omparable to H 2 cooling , 
as long as 10 3 K < T < 10 4 K dSusa fc Umemurd[200l . 
Therefore, radiative cooling by heavy elements in these 
halos do not play central role as long as complete mixing 
is assumed. 

On the other hand, metals in the interior of expanding 
SNR will mix with the shocked shell it sweeps up be- 
cause of Rayleigh-Taylor and Kelvin-Helmholtz instabil- 
ities. This will enrich the shocked region to much higher 
metallicities than Z/Zq ~ 10 -35 , radiatively cool the 
shell much faster than any H2 it has formed or swept 
up, and then fragment clumps. It may be a clump that 
collides with the neighbor halo, not an intact shell, but 
this is for future numerical studies to determine. 

3. EFFECTS OF SUPERNOVA SHOCK ON THE NEIGHBORING 
DARK HALOS 

In this section, we assess the effect of the SN shock onto 
neighboring dark halos by analytic arguments comparing 
timescales. 

3.1. Gas evacuation by shock momentum 

In case a neighboring halo is located very close to the 
center of the SN explosion, the gas in the halo would be 
evacuated by the shock momentum. As a result, subse- 
quent star formation in the halo should be inhibited. As- 
suming the momentum conservation after the gas behind 
the shock cools, we can roughly evaluate the velocity of 
the gas in the halo after the impact of the shocked shell. 

In Figure [TJ we represent three cases of gas evacua- 
tion by shock momentum. First, we consider the simple 




Fig. 1. — Three cases of gas evacuation by shock momentum 
are represented in this figure. The shock propagates at tig before 
the shocked gas cools by radiative cooling. After that, the shock 
propagates under the momentum conservation. The intermediate 
case (case(ii)) indicates that the shock shifts to the momentum- 
driven expansion phase in the halo. 



momentum conserving case (case(i) in Figure [lj, where 
/ CO oi < 2r h (i.e., t cool < 2r h /t>s) is satisfied. Here, / coo i 
denotes the distance that the shock sweeps the halo be- 
fore it enters the momentum conserving phase, i.e., 

/cool = min (v s tcool (Tps), 2r h ) . (11) 

Here, is the radius of the dark halo. The equation 
of the momentum conservation is given as follows: 

( n b \ 

m s v s = m s + 7— M dh Vtot, (12) 
\ "M / 

where Mdh is the mass of the dark halo and mg is the 
mass of the gas shell colliding with the halo, which is 
given as 

m s = (Mpopm + ^R s 3 Po^j (13) 

Here, Mpopm is the mass of progenitor POPIII star 
and we assume that the whole mass of this POPIII star 
is released as ejecta. We assume Mpopm = 14OM in 
the present calculation. 

On the other hand, in case of t coo \ < 2r^/vs (i-e., 
/cool < 2rh, case(ii) in Figure [1]), the gas in the halo 
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Fig. 2. — Velocity of the shocked gas (vtat) is plotted as a function 
of the distance from the SN center, ilg, for i?SN = 10 52 erg and 
M dh = 2 X 1O 5 M0. The dotted line represents the escape velocity 
of the halo with M dh = 2 X 1O 5 M . 



Fig. 3. — Condition of gas evacuation by the SN shock for 
E>SN = 10 52 erg. The horizontal axis shows the distance between 
the nearby halo and SN center, while the vertical axis denotes the 
mass of the halo. In the hatched area («tot < Vise), the shocked 
gas is bounded inside the dark halo. The unhatched region labeled 
as vtot > Vcsc, the gas will be evacuated by the shock from the SN. 



is cooled while the shock propagates in the halo by ra- 
diative cooling, and it evolves into the momentum con- 
serving phase. Since the gas pressure of the halo be- 
hind the shock front can keep pushing the shock before 
it cools, Vtot cannot be estimated by simple momentum 
conservation equation. It is difficult to assess this effect 
analytically. Thus, we assume that the shock velocity 
is not affected by the halo while the radiative cooling is 
still inefficient even after the shock enters the halo. This 
assumption could be oversimplification, since even adia- 
batic shock will slow down to some extent when it enters 
the dense region. In order to quantify this effect, numer- 
ical study would be necessary, which is beyond the scope 
of present study. Therefore, we have to keep in mind that 
the shock heating effect is maximally taken into consid- 
eration in this model. In this intermediate case ((ii) in 
Figure [1} , we assume that the momentum conservation 
after shocked gas in the halo is cooled: 

(ms + ™° l b M dh ) i>s = (m s + J «tot(14) 

V 2r h S2 M / \ "m / 

In case the shock propagate through the entire halo 
within the cooling time of the shocked gas, that is the 
case of /cool = 2rh, the shock velocity in the halo is as 
large as vs during the shock propagation because of the 
inefficient cooling (case (iii) in Figure [TJ). 

The velocity of the shocked gas shell, Vtot, is obtained 
from the equation (TT4|) as follows: 

. ^COOl ii j 

«tot = o vs- (15) 

m s + 7=j— Mdh 

This equation includes all of the three cases. The cases 
of icooi ^ 2rh/ws and / coo i = 2r^ are the limits of effi- 
cient/inefficient cooling. 

The gas will be still bounded in the halo even after 
the shock arrival, in case vtot is smaller than the escape 
velocity V csc of the halo. On the other hand, if u tot > 
V csc , the gas is blown away from the halo potential, in 
which subsequent star formation in the halo is totally 
quenched. Thus, we regard the condition v to t < V csc as 



a necessary condition for the triggered star formation in 
the neighboring dark halos. 

In Figure [5J Vtot , as a function of the distance from 
the SN center assuming Mdh = 2 x 10 5 M Q and E$n — 
10 52 erg. In this case, Vtot is smaller than V esc in the 
ranges 0.7kpc < R s < l.lkpc and R s > 2.3kpc. There- 
fore, in case the neighboring halo is located at such dis- 
tance, the gas in the halo is not evacuated by the shock 
momentum. 

Figure shows the gas evacuation from the halos for 
Esn = 10 52 erg. The vertical axis denotes the mass of 
the dark matter halos (Mdh), whereas the horizontal 
axis represents the distance from the SN center. In the 
hatched area denoted as v to t < V csc , the gas in the dark 
halos are not lost by the shock momentum. It is worth 
noting that the hatched area around R$ = lkpc corre- 
sponds to Tps ~ 10 4 K, where the Lya cooling dominates 
the others. 

3.2. Cooling/collapse of the shocked gas 

The halos that satisfies the condition i> to t < V csc are 
able to survive the disruption by the SN shock momen- 
tum. As a next step, we consider whether those survived 
halos can collapse or not. We set the collapse criteria 
as the condition where the gas temperature is decreased 
below T v i r , before the gas cloud expands by the thermal 
pressure, or bounces by adiabatic compression. In other 
words, if following conditions 

t coo i{T) < t s (l coo i) and t cool (T) < t s 

are both satisfied until they cool below the initial virial 
temperature, we regard that the gas in such halos can col- 
lapse to form stars. Here, we consider the sound cross- 
ing timescale, t s (i C ool) — ^cooi/cs, to be the expansion 
timescale of the shocked region. In general, the con- 
dition t coo \ < iff is regarded as the collapse criterion 
of gr avitationally bound system (e.g., iRees fc Ostrikerl 
1977). Once this condition is satisfied, the system starts 
to collapse. The cooling timescale in low-density gas is 
basically inversely proportional to the density, while the 
free-fall timescale proportional to the inverse of square 
root of density. As a result, the ratio icooiAff gets smaller 
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Fig. 4. — Evolution of timescales in the shock-heated gas is shown. Four panels correspond to various distances Rg (upper-right corner), 
where the dark halos are located. t coo i, t B , iff represent the cooling timescale (solid), the sound crossing timescale of the shock- heated gas 
(dashed) , and the free-fall timescale (dotted), respectively. The thick dot-dashed line indicates the virial temperature of the dark halo in 
the case of M va ^ = 2 X 10 5 Mq, whereas the thin dot-dashed line shows the postshock temperature Tpg. 



as the collapse proceeds. That 's why the condition is re- 
garded as the collapse criterion. In case we consider the 
primordial star formation, the cooling timescale in most 
of the final run-away collapse phase is not proportional 
to the inverse of density, however, the absolute value of 
the cooling timescale is shorter than the other timescales 
outside the halo. Thus, the collapse of the cloud contin- 
ues following the track along which t coo i = tg is satisfied. 
In addition, we have to consider the adiabatic expansion 
of the shocked gas, since the sound crossing timescale of 
the shocked gas could be very short in the present case. 
Thus, we need to add the inequality t coo \ < t s (l coo \) to 
the condition of the collapse criteria. 

In order to assess this condition, we have to follow the 
thermal evolution of the shock heated gas in the dark 
halo. Each of the four panels in Figure [4] shows the evo- 
lution of timescales after the shock heating for halos with 
M dh = 2 x 10 5 M Q located at R s = 250pc, 480pc, lOOOpc, 
and 1300pc. Thermal energy of the SN explosion is 
assumed as -Esn = 10 52 erg. The horizontal axes de- 
note the temperature of the gas on the way of cool- 
ing, whereas the vertical axes show timescales. In or- 
der to assess the cooling timescale for T < 10 4 K, we 
have to take into account the H2 cooling, which is pro- 
portional to H2 fraction. H2 fraction could be obtained 
by solving nonequilibrium chemical reaction equations, 
however, approximate values are already obtained. We 
use /h 2 = 10~ 4 in case the postshock t emperature, Tps, 
is less than 10 4 K (jNishi fc Susalll999D . whereas /h 2 = 



2 x 10~ 3 is employed for Tpr > 10 4 K (lOh fc Haimanl 
[200l iShapiro fc Kandfl987l : iSusa et aTlll998h . 

If the halo with M d h = 2 x 1O 5 M is located at 
R s = 250pc (upper- left panel), the shocked tempera- 
ture is too high for the gas to remain the gas inside the 
halo potential. In fact, expansion timescale is already 
shorter than the cooling timescale at T — Tpg. Thus, 
the gas in this halo is lost because of the shock heat- 
ing. On the other hand, in the cases of Rs — 480pc 
(upper- right panel), Rs = lOOOpc (lower-left panel), and 
Rs = 1300pc (lower-right panel), the gas in the halo can 
start to cool because i C ooi(T) is smaller than t s Gcooi) at 
T = Tps- Among these examples, t coo i(T) is shorter 
than ts and t s (l C ooi) as long as T > T V1Y is satisfied in 
the case Rs = 480pc. Therefore, the gas can be cooled 
below T v j r before it expands, and the cooling phase could 
be followed by gravitational contraction. In contrast, in 
the cases of Rs = lOOOpc and Rs = 1300pc, the cooling 
process becomes inefficient before T = T vlr is achieved. 
In the case of Rs = lOOOpc, the postshock temperature 
exceeds 10 4 K. Thus, /h, = 2 x 10~ 3 is achieved even for 
T < 10 4 K, whereas / Ha = 10~ 4 for R s = 1300pc. There- 
fore, the cooling condition is not simply determined by 
the efficiency of H2 formation, although it is necessary 
for the cooling below 10 4 K. 

The allowed region for cooling/collapse of the shocked 
gas is shown in the next section. 



7 




12 3 

Rs [kpc] 



Fig. 5. — Regions in which the shocked gas is cooled and the 
gas in the dark halos is not evacuated on the Rg — M^h plane. 
Three panels denoted as (a), (b), (c) represent Esn = 10 51 ,10 52 , 
and 10 53 erg, respectively. In the hatched area labeled as "Survive" 
(i'tot < Vesc), the shocked gas is bounded in the halo potential. On 
the other hand, in the shaded region denoted as "Cool" , postshock 
gas cools below T v ; r . 



4. POSSIBILITY OF TRIGGERED STAR FORMATION IN THE 
NEIGHBORING DARK HALOS 

Now, we are ready to combine all conditions to find 
the criteria for positive/negative feedback effects by first 
SNe on the neighboring dark halos. Three panels in Fig- 
ure \5\ correspond to the results with Esn = 10 51 ,10 52 , 
and 10 53 erg, respectively. Extremely energetic cases with 
10 52 and 10 53 erg can be interpreted as the hypernova or 
the pair-instability SN. The hatched area bounded by 
dotted curves (labeled as "Survive" ) denotes the regions 
in which i> to t < V csc is satisfied. The shaded region de- 
noted as "Cool" represents the condition in which the 
cooling rate is high enough for the gas in the halo to 
proceed further gravitational contraction. 

First of all, gas components in dark matter halos con- 
sidered here (M ^ 10 7 Mq) are blown away by the shock 
momentum if the halos are close enough to the SN center. 
The critical distance within which the gas is evacuated 
falls within the range ~ 0.3 — 1.5kpc depending on the SN 
energy and mass of the halo. In case we assume normal 



core-collapse SN (E'sn = 10 51 erg), the critical distance 
is ~ 0.3 — 0.5 kpc, whereas it is ~ 1 — 2kpc for pair- 
instability SN (Esn = 10 53 erg). The mass dependence is 
not so strong, but basically the low-mass halos are more 
fragile than the massive halos. In addition, we also note 
that the boundaries of the 'Survive' regions have com- 
plicated structures at Rs > lkpc, reflecting the shape of 
the cooling function. Recent studies suggested that halos 
further than ~ 1.5kpc from the original star will not be 
reached by the blast. In the present paper, we include 
distances greater than this in our analysis, but we are 
doing so only for completeness. 

The region where "Survive" and "Cool" are compatible 
with each other correspond to the case that the shock- 
heated gas successfully cooled and collapse in the dark 
halo. However, it is clear that the halos which can sur- 
vive the evacuation by the SN shock and collapse do not 
exist. Consequently, the SN feedback has basically neg- 
ative effects on the star formation in surrounding halos. 

5. DISCUSSION 

The nearby dark halos come under the influence of 
the radiation from the SN progenitor star on the satel- 
lite halos prior to the explosion. Thus, many low-mass 
dark halos might be photoevaporated by the UV flux 
of progenitor star. One-dimensional r adiation hydro- 
dynam ics simulations are performed by Aim & Shapiro 
(|2007f ) on this issue. They found that effects of radiative 
feedback on the gravitational contraction of the gas in 
low-mass halos (Mdh < 1O 6 M0) are very complicated. 
The feedback effects are qualitatively different depend- 
ing on the distance from the source star, evolutionary 
stage of the halo at the onset of r a diative feedback, and 
mass of the halo. iWhalen et al.l (|2008af) also have in- 
vestigated on this problem by two-dimensional radiation 
hydrodynamics simulations, for a single halo with mass 
of Mdh = 1-35 x 1O 5 M , in which they take into account 
a detailed evolutionary stage of the halo. They found 
that most of the gas in this halo is photoevaporated by 
the SN progenitor star prior to its death at very early 
evolutionary stages of the halo. On the other hand, they 
indicated that the /-front could not reach the core of the 
halo and this core survived quite well if central d e nsities 
of the halo rose beyond 50cm~ 3 . lYoshida et al.l (|2007l ) 
also performed three-dimensional radiation hydrodynam- 
ics simulations with realistic cosmological density field. 
They found the H n region of 100M o POPIII star extend 
to ~ lkpc, within which most of low-mass halos are pho- 
toevaporated. However, the number of survived halos 
through UV flux from progenitor star is still under de- 
bate, since we only have the results by multidimensional 
simulations with restricted parameter space (mass of the 
halo, mass of the progenitor star). In any case, it is worth 
to investigate the effects of SN shock on the nearby halos 
with various parameters by analytical calculations. 

In our study, we completely ignore mass loss in the 
neighbor halos due to photoevaporation by the progen- 
itor star. Therefore, it must be noted that our results 
are taken to be a lower limit on the damage done by 
the expanding remnant to the halo because the halos the 
SN blast actually encounters will have less gas. Also the 
negative feedback effect we have found applies only to 
the one scenario, in which the SNR increasingly sweeps 
up mass and then collides with the neighboring halo. It 
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could be that SN directly form more stars by other means 
than those they quench in nearby halos in our mecha- 
nism. 

6. CONCLUSION 

We have studied the feedback effects by first SNe with 
Esn = 10 51 , 10 52 , and 10 53 erg on their neighboring dark 
matter halos. Consequently, the conclusion can be sum- 
marized as follows. We find that the star formation activ- 
ities in the neighboring dark matter halos (M 10 7 M©) 
are basically suppressed in case they are located close to 
the SN center, because of the gas evacuation effect. The 
critical distance within which the gas is blown away falls 
within the range ~ 0.3 — 1.5kpc depending on the SN 
energy and the halo mass. In case we assume normal 
core-collapse SN (-Esn = 10 51 erg), the critical distance 



is ~ 0.3 — 0.5 kpc, whereas it is ~ 1 — 1.5kpc for pair- 
instability SN (Esn = 10 53 erg). In addition, we find 
there is very little window in the parameter space where 
star formation activities in dark halos are induced or pro- 
moted by neighboring SN. 
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